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Abstract 

The PI control of first-order linear passive systems through a delayed 
communication channel is revisited in light of the relative stability con¬ 
cept called (T-stability. Treating the delayed communication channel as 
a transport PDE, the passivity of the overall control-loop is guaranteed, 
resulting in a closed-loop system of neutral nature. Spectral methods are 
then applied to the system to obtain a complete stability map. In particu¬ 
lar, we perform the H-subdivision method to declare the exact (T-stability 
regions in the space of PI parameters. This framework is then utilized 
to analytically determine the maximum achievable exponential decay rate 
(t 2 of the system while achieving the PI tuning as explicit function of (tJ 
and system parameters. 


1 Introduction 


Passivity-based control relies on the fact that the power-preser ving intercon¬ 


necti on of two passive subsystems yields again a passive system [Ortega et al 
I998l| . Lyapunov stability of the interconnected system follows from passiv¬ 


ity while asymptotic stability is usually achieved by adding appropriate damp¬ 
ing. Due to its simplicity and robustness, passivity-based control has attracted 
rese archers and practitioners in the c on trol community for several decades 
e.g. lYoula et ah . I959l WillemsI . 1972 al ibi. Hill and Movlan . 1976 . Byrnes et al 


19911 


van der Schaftl . 200Cll| . 


However, if a delayed communication channel stands between the plant and 
the controller, as in a typical control scenario, passivity arguments fail due to 
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non-passive properties of the channel. Here, the loss of passivity follows from 
the fact that the Nyquist plot of a pure delay does not lie in the right-hand 
side of the comple x plane. In their seminal w ork, inspired from the study of 
transmission lines Anderson and Soond 1989l| proposed a useful modification 
of the communication channel to remedy the aforementioned design problem. 
More precisely, the communication channel is transformed into a passive system, 
thus recovering the simplicity and effectiveness of the passivity-based design. 
This idea has been discussed in many contributions and has given rise to an 
outstanding number of proposals addressing optimality issues, applications in 
th e field of robotics, motor control, among other studies. The reader is referred 
to iNuno et al.l 2011 1 for a recent tutorial. 

This paper revisits the modified communication channel from the perspec¬ 
tive of time-delay systems theory using spectral methods considering first-order 
linear passive systems. The problem is motivated by the wide variety of indus¬ 
trial processes described by first-order plants with time-delay and commonly 
regulated by PI controllers Silva et al.l . 12001 1. such as DC servomotors exten¬ 
sively used in industry. In this paper, the emphasis is put on the performance 
of the closed-loop system when a communication channel stands between the 
plant and the controller, which is quantified by its cr-stability degree. Here, a 
approximates the exponential decay rate of the system response. 

Our analysis is b ased on classical resu l ts of t ime-delay systems of r etarded 


and neutral nature [Bellman and Cookd . 119631 iHale and Lunell . Il993j . The 


problem under consideration, though infinite dimensional, involves a reduced 
number of parameters. Hence, a comprehensive frequency domain analysis of 
the closed-loop characteristic quasipolynomial can be perform ed. Parti c ularl'v 


rtic uiariy , 

3[i9i3, 


we depl oy a critical extens ion of the ^-subdivision method of iNeimarl ^ 
see also Sioahi et al.l 20111 for advanced methods, which consists on (i) the de¬ 
termination of the stability boundaries corresponding to roots at s = —a and 
s = —a + ju}, which provides a partition of the space of PI parameters and (ii) 
the verification of the relative stability degree a of each region in the partition. 
Having generated the complete set of cr-stability boundaries and determined 
the cr-stability regions, the exact cr-stability maps follow. This framework then 
results in a fully analytic characterization of the maximum achievable expo¬ 
nential decay rate for which simple tuning formulae for practitioners are finally 
declared. 

The contribution is organized as follows: In Section |2l the delay-free and 
the fixed, non-zero delay cases are analyzed, illustrating the failure and loss of 
performance in the passivity-based design strategy. In Section |31 the scattering 
transformation is introduced and the characteristic function is obtained, the cr- 
stability maps are sketched and tuning rules for points of interest are derived. A 
theoretical limit on the cr-stability is found and a tuning rule for the scattering 
transformation is proposed. It is shown that, when the rule is followed, the 
theoretical limit can be approached arbitrarily close. Concluding remarks are 
given in Section |4l 
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Figure 1: A control-loop including a delayed communication channel. 


2 Problem statement 

Consider a first-order linear system of the form 


X = —ax + bui , (la) 

yi=x , (lb) 

where iti, j/i and a; £ M are the input, output and state, respectively. The 
parameters a and b are assumed to be non negative, which ensures passivity 
with storage function Vi{x) = x^/{2b). 

Consider the PI controller 

i = uo , (2a) 

2/0 = kpUQ + ki^ , (2b) 

where uq, yo, ^ € M. are the controller input, output and state, respectively. The 
proportional and integral gains kp and ki are assumed to be positive, so the 
controller is also passive. The system and the controller are interconnected as 
per the following pattern 


ui = - 2/0 and uq = yi - y/. 


Since the interconnection of two passive systems is again passive, the closed- 
loop characteristic polynomial is stable for all gains. 

However, when a communication channel with delays is introduced in the 
loop, as shown in Fig. [1] the closed-loop transfer function takes the form 

2 / 1 ( 5 ) ^ {kpS + 

yl{s) + as + {kpS + ki)be~^^ 


with hi the forward delay, /12 the return delay and h = hi + h 2 the round- 
trip delay. The stability properties of the closed-loop are then defined by the 
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location of the roots of the characteristic equation 

p{s) = + as + {kpS + ki)be~’^^ , (4) 

also known as the characteristic quasipolynomial. Notice that the presence of 
the delay h > 0 in the communication channel induces infinite-dimensionality 
to the system due to the exponential term and therefore the quasipolynomial 
in (jll) bears an infinite number of roots. Since it is impossible to compute all 
these roots, the stabilization of the zero-solution is not trivial. Moreover, since 
the delay channel is not passive, the passivity argument fails and stability can 
no longer be ensured for every combination of positive parameters. 

In the following we consider the problem of finding the setting on the pa¬ 
rameters (fcp, ki) that create the maximum decay rate for the system (II|)-(I11) in 
the presence/absence of time-delays. 


2.1 Performance degradation as a result of delays 


As a preliminary step for the characterization of the maximum decay rate, we 
begin with the decomposition of the (fcp, fci)-plane. Besides pure stability, we will 
be concerned with the exponential decay of system solutions with a given degree 
cr, that is, with the cr-stability of the system. This happens only if all the roots o f 
the characteristic quasipolynomial have real parts less than — cr Gu et alll2003l| . 
Then, we will determine the set of all {kp,ki) points for which the closed-loop 
system is cr-stable. 

Following the ^-subdivision method, we equate (|4]) to zero and set s = —a, 
which gives the boundary 


h = ak„ + 


crag 

bf^ho 


(5) 


Now, setting s = —a + juj and solving for kp and ki gives the parametric 
equations 


kp{uj) 


ki{uj) 


—a 2 aOJ cos{hoj) -I- [aa„ + sin(hw) 

(cr^ -|- (w cos{huj) + Qa sm{huj)) 
buje^'^ 


(6a) 

(6b) 


Equipped with ([5]) and under continuity arguments, we can now declare the 
exact regions in (kp, ki)-plaiie for which the quasipolynomial Q is cr-stable. 

Remark on non-dimensionalization. In investigating the spectral properties of 
the considered control-loop it is not necessary to distinguish between every 
possible combination of (a, b) parameters. In fact, introducing the following 
scaling variables t = at and u'l = (b/a)ui, the non-dimensional form of the 
general plant m is obtained as 


^x(t) = -x{t) +u[{t) , 
yi ( r ) = x{t) , 
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where the dependence on (a, b) is obviated. In other words, without loss of 
generality we can set a = & = 1 in the figures that follow for demonstration of 
the control-loop properties. In consequence, and consistent with the above non- 
dimensionalization, the figures for all possible (a, b) combinations might change 
quantitatively, but the qualitative features remain unchanged. 


Considering that a is given, equations ©, © decompose the {kp, fci)-plane 
into a finite number of disjoint regions. Due to continuity arguments, each of 
these regions is then characterized by the same number of strictly cr-unstable 
characteristic roots. We will refer to the collection of all pairs {kp, ki) for which 
the number of cr-unstable roots of p{s) is zero as the cr-stability domain T)^. 

Using a = b = \ in ©, ©, we obtain Fig. [2]for h > 0. In order to assist 
the reader, for a given a the corresponding Va is filled in color and the stability 
crossing boundaries are trimmed to fit the boundary of 'D^. Two interesting 
observations are in order. Firstly, in Fig. |2(b)| one can observe that searching 
for a faster response will eventually result in the collapse of T>^. We have that, 
at the critical value a* where the boundaries obtained from © , © meet each 
other at Do-*, a triple real dominant root will be generated at —a*. This point 
corresponds to the maximum achievable closed-loop exponential decay rate cr*, 
and these argume nts are consistent with previou s discussions on the maximum 


achievable a, see [Michiels and Niculescul . l2007j| Theorem 7.6. Secondly, from 


Fig. |2(a)[ in contrast with the delayed case, it is clear that achieving an arbi¬ 
trarily large exponential decay a is possible in the absence of communication 
delays. 

Let us characterize the maximal achievable decay cr* for the delayed case. 
This characterization is then particularized to conclude on the relative stability 
properties of the delay-free case. 


Proposition 1. Consider a plant o in closed-loop with a PI controller © 
satisfying kp > 0 and ki > 0. A delayed communication channel with round-trip 
delay h > 0 stands between the system and the controller. 


(i) The maximal achievable exponential decay is given by 


a 


* 


4:-\- ah — \/8 + 
2h 


(7) 


(ii) The minimal PI controller gains assigning a given cr* > a > aj^ are 


kp — 


ki = 


ahOrr — 02(7 
bf,ha ■ 

(T'^{hac -\- 1 ) 

bf,hcr 


(8a) 

(8b) 


where condition a > a/2 ensures that kp > 0 for all h > 0. 

Proof. (i) For a given a, the corresponding stability domain is delimited by 
the parametric equations corresponding to a real root at s = — tr and a 
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(a) Delay-free case, h = 0. Achievable decay is unbounded. 



Figure 2: cr-stability regions. Communication channel without scattering trans¬ 
formation. The diamond markers correspond to minimal gains ensuring a given 
a. 
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pair of complex roots s = —a ± juj of the quasipolynomial (|3]). As a 
consequence, the collapse of the cr-stability domain occurs at a triple root 
at —cr. Thus, the quasipolynomial dl]), its first and its second derivatives 
must vanish at s = —a. That is, 

—bae^'^kp + be^'^ki = aa^ , (9a) 

+ ha)kp — hbe^'^ki = —a^a , (9b) 

—hbe^’^{2 + ha)kp + h^be^'^ki = —2 . (9c) 

Solving these equations for a gives 

4 + aft. ± Vs + a^ft^ 

=-2ft- ■ 

The solution that ensures kp > 0 and fti > 0 is ©• 

(ii) The minimal gains occur at the intersection of the two boundaries. Hence, 
the conditions (15a)) and (I9bp must hold. The result follows by solving these 
equations for kp and ki. 

□ 

We end this section with a comment on the performance degradation as a 
result of delays. To this end, let ft —>■ 0 in ((T]). It follows that the maximum 
achievable decay rate tends to infinite as the delay vanishes. From a practical 
point of view, an upper bound for a* is solely determined by the physical limi¬ 
tations of the considered delay-free system. Moreover the minimal PI controller 
gains assigning a given u > a/2 are given by kp in (l8aD and by ki in (I8bl) with 
ft = 0 . 

Finally the rest of the paper is dedicated to the problems of (i) improving 
the relative stability of the system when challenged by the presence of commu¬ 
nication delays (ii) recovering the passive properties of the overall control-loop 
and (hi) algebraically designing the controller gains to prescribe a desired ex¬ 
ponential decay rate. 


3 The scattering transformation 


A classical approach to the study of transmi ssion lines co nsists in applying a 
linear transformation on the state variables Chenel . Il992 |. By applying such 
transformation, the transmission line equations, i.e., the telegrapher’s equations, 
transform into a pair of uncoupled delay equations, i.e., transport PDF. It is 
then possible to understand the dynamics of the transmission line in terms of 
wave propagation. 

The reverse argument was proposed by Anderson and Soonel 1989l | : Suppose 
we have a communication channel consisting of a pair of delays. Apply the 
inverse transformation to emulate the behaviour of a transmission line. Since 
transmission lines are passive (lossless), the passivity argument is restored. 
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More prec isely, consider a pair of delays given by the transport PDE, 
Krstic . 2009l| for applications of transport PDE to backstopping design, 


see 


fs+{l,t)\ 0\fs+il,t)\ 

\Si {l,t)J 0 h2j \st il,t)J 


( 10 ) 


where are the partial derivatives of the scattering variables with respect 
to the spatial variable I G [0,1]. Notice that, at the boundaries, the solu¬ 
tions satisfy s+(l,t) = s+(0,t — hi) and = s~{0,t + /i 2 ). This is the 

communication channel. Consider now the linear transformation 


^ A d\ ^ fs+{l,t)\ 
\v{l,t)J 1^1 -dj \s-{l,t)J 


( 11 ) 


where d > 0 is a design parameter. It follows from m and (US, with h — 
hi + h 2 , in general hi ^ / 12 , that 


(^ _h ( 0 d\ / 

2 \l/d Oj ’ 


where and (/it,r't) are respectively the spatial and temporal partial 

derivatives of {fi, v). The above PDE corresponds to the (lossless) telegrapher’s 
equations. 

Transformation (ED can be enforced at the boundaries of ED- However, 
it is necessary to be cautious in respecting the causality of the system: The 
variables s(j" := s“''(0,t), s(" := s~{l,t), /tq := /r(Oji) '*-'1 •= '*-'(1)^) are free, 
while := s'''(l,t), Sq := vq := v{0,t) and fii := fii{l,t) depend on 

the former variables and their past values. With these restrictions in mind, we 
can set, at the boundary I = 0, 



(cf. (ITT]) at ^ = 0). At Z = 1 we have 



-{“) (A ■ 

!) (4) ■ 


(cf. (fTTl) at Z = 1). Here, equatio n (IT^ is known as the scattering transformation . 


Further details can be found in Nufio et ah . 2011 . 20091 Niemever and Slotine . 


|2004 [l9^ |. Finally, we can use the interconnection pattern 


as in FiglS] 


uo = vo-yi , yo = -yo , ui = fj-i and vi = yi , 


(13) 


















PI controller 


Passive system 




iti 


Modified communication channel 



yi 


Figure 3: A control-loop including a communication channel with the scattering 
transformation. 


3.1 Recovering stability and improving a-stability 

When the scattering transformations are introduced as shown in Fig.[3J and the 
parameter d is arbitrary, the closed-loop transfer function is obtained from o, 
m and m and using both the interconnection pattern (IT^ and the equations 
of the delayed communication channel = SQ{t — hi) and s,J = Si{t — h 2 ). 
Then, we have 

yi(s) 2d{kpS + ki)be~^^^ 
yt{s) P 2 {s)s'^ + pi{s)s + pq{s) ’ 

where 

P2{s) = {l + e-’^^)d+{l-e-'^^)kp, 

pi{s) = (1 -b e~^‘‘){bkp + a)d -b (1 - -b akp -b h) , 

Po(s) = (1 + e~'^‘‘)bkid -b (1 - e~^'‘)aki . 


The characteristic quasipolynomial is thus 

p(s) =P2(s)s^+Pi(s)s+Po(s) • (14) 


Notice that this quasipolynomial is of neutral type. In other words, its pr i ncipal 


coefficient, P 2 (s), contains ex ponential terms, see Bellman and Coom Il963l 


Kharitonov and Mondi4l2011 | for a discussion and analysis of quasipolynomials. 


In the time domain, this means that the time derivative of the state of the system 
does not only depend on delayed states, but also on the derivative of the delayed 
state. 

To determine the boundaries of the cr-stability region, we equate the quasipoly¬ 
nomial m to zero, set s = —a and solve for ki as a function of kp, 


{1 + + {1 - e^^)bd 

* ^ ^ ^ (1 _ e^®)acr + (1 + e^’^)bd 


(15) 
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Now, setting s = —a+juj and solving for kp and ki gives the implicit parametric 
equations 

^(w) (^^) = d • S(a;) , (16) 

where 

Aii{uj) = (qw — aa)bd — f 3 {aaa + uj^) — "fuja2cr , 

^12 (w) = abd + Paa — jw , 

^21 (w) = (7CT + aoj)bd — ^{aacy + + Puja 2 a , 

^22 (w) = -^bd + qocr + Poj , 

Bi{uj) = {Pa + ')Oj)bd + a(aaa- + uj^) — 'yuja 2 cr , 

B2{to) = (qcr — Puj)bd — 'y{aacr + — aijja2a ■ 

and 

a = 1 + 6^“^ cos{huj) , 

P = I — cos{huj) , 

7 = sm{huj) . 


Recall that a necessary condition for the stability of neutral-type delay sys- 


terns is the stability of the difference operator iHale and Lunell . Il993liMichiels and Niculescu 


2007 . Olgac and Sioahi . 2004 , Olgac et ah . 2008l| . which is the inverse Laplace 


transform of P 2 {s)- When tr = 0, the difference equation becomes 


{d + kp) x{t) + (d — kp) x{t — h) = Q . 

Since \{d — kp) / {d + kp)\ < 1, the stability of the difference operator always 
holds. On the contrary, when ct > 0, the characteristic equation of the difference 
operator becomes d + kp + e“^®e^'^(d — kp) = 0. The necessity on the stability 
of the difference operator imposes the new condition 


e 


h(7 


d — kp 
d -\- kp 


< 1 


or, equivalently. 




- 1 


_/l(7 


< k„ < d- 


+ 1 


(.ha ^ '-P ^ “gfta _ I ■ 

The expressions (HSl), dm) and (II3 are used to determine the cr-stability re¬ 
gions in the {kp, /ci, d)-space of parameters. These regions are shown in Fig. |4(a)| 
For illustration purposes, the slice d = 15 is also shown in Fig. |4(b)| 


Remark 1. The scattering transformation recovers the striking property, ob¬ 
served in the delay-free case, that the whole first quadrant of the parameter 
space ensures stability of the closed-loop. Regarding cr-stability, the size of the 
regions in Fig. |4(b)] are larger than those in Fig. |2(b)] Also, the maximal achiev¬ 
able exponential decay, a*, is greater with the scattering transformation than 
without it for sufficiently large d. 
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(a) Regions of (j-stability in the (/cp, fcj, d)-space of parameters. The upper 
parts of the regions have been clipped so that inner regions corresponding to 
higher cr are visible. The planes d = kp and d = Cmin^p shown in light 
gray. 



(b) (T-stability regions on the slice d = 15. The diamond markers correspond 
to the minimal gains that ensure a given cr. Maximal achievable decay is 


Figure 4: Communication channel with delays, h = 0.1, and scattering trans¬ 
formation. 
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Again, we characterize the maximal achievable decay and the corresponding 
control gains. 

Proposition 2. Consider a plant CD in closed-loop with a PI controller m 
satisfying kp > 0 and ki > 0. A communication channel with round-trip delay 
h > 0 stands between the system and the controller. For a given d > 0, apply 
the scattering transformations m at the channel end points. 

(i) The maximal achievable decay is a root a of md{<j), where 


md{cr) = (l+e'“") {2hb^d^ + + 4)b'^d'^ + 2h{a^ - al^)bd - h^aal) 

+ (1 — {hfab^d^ + 2h{a + a)b‘^d^ + (4a — h'^aa'^)bd — 2ha^) (18) 


such that the gains 


{bd - + 2ae'“" (2bd + hlpd"^ - al)) - (bd + 

kp = d - p - 

{bd + 0(7 + e^'^{bd — Oa)) 

k, = ^ 

{bd + Ocr + e^'^fbd — Oa)) 


, (19a) 
(19b) 


are non negative and such that C3 holds. 


(ii) The minimal PI controller gains assigning a given a > al2 are given 

by CID- 

Proof. (i) As in the case where no scattering transformation is employed, the 
cr-stability regions collapse at a triple root at s = —a. The fact that the 
quasipolynomial (jl4[) and its first derivative are zero give the conditions, 
written in compact form. 



where 

All = —(1 + e^'^)abd — (1 — e^'^)aaa , 
yli2 = (l + e'‘")M+(l-e'‘‘^)a<7 , 

A21 = bd + 02(7 + 6^’’^ {{^d — acr){h(T + 1) + cr) , 
A22 = 1 - e^'^{h{bd — Ocr) + 1) , 


and 


i?i = (1 — e^'^)abd + (1 + e^'^)aacr , 

i?2 = -{bd-\- a2cr) + {{bd - ac7){ha + 1) + ct) . 
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The solution of this linear system of equations is given by (HU). Computing 
the second derivative of da gives 


(A31 A32) =d-i 33 , 

where 


( 21 ) 


A 31 = (2 - {{h{bd - ttcr) + l){ha + 2) + ha)) , 

A32 = he^^{h{bd-aa)+2) , 

B 3 = - (e^'^ {{h{bd - a^) + l){ha + 2 ) + ha) + 2 ) . 

Substituting (ITOl) into m gives the implicit equation md{a) = 0. Then, 
is given as a root of md(a) satisfying the non negative condition on the 
PI control gains and the necessary stability condition da associated to 
the neutral nature of the quasipolynomial. 

(ii) The minimal control gains occur at the intersection of two boundaries. 
The result follows by noting that at this intersection a double root of the 
quasipolynomial da arises, in other words, when (IT^ holds. 

□ 

Observe that the vanishing of the quasipolynomial (fl3 and its first two 
derivatives are necessary conditions for the tr-stability regions to collapse. There¬ 
fore, the roots of da must be verified via back substitution into the solution 
of equation (1201) . Then, if fcp > 0, > 0 and da hold, the CT^-stability of the 

detected collapse point has to be checked for tr^ to be feasible. 

Remark 2. It follows from the implicit function theorem that crj can also be 
characterized as the value of a at which the derivative of (|19ap with respect to 
a is equal to zero. 

3.2 Least upper bound on the exponential decay rate 

To provide an idea of how restrictive (I17p is, we have computed the minimal 
gains (fT^ and plotted them in Fig. [5] Notice that, when the minimal gains are 
used, the upper bound in (1171) is never infringed. The restriction a > a/2, on 
the other hand, ensures the lower bound. 

It is clear from Fig. [5] that, as d —>■ 00 , the maximal exponential decays 
a/i accumulate at a point, which we denote by (Tsup- This is formalized in the 
following proposition. 

Proposition 3. The maximal exponential decays a^ are bounded by CTsup; which 
is defined implicitly by 

2(1 -I- -I- /icrsup(l - = 0 . (22) 
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Fiffiire 5: Minimal k„ as a function of a (llQaD. normalized bv d. The bounds 03 
are also shown. The diamond markers correspond to cr^, in accordance with (jl8l) . 
Notice that kp{a'^)/d converges to a constant value as d —>■ oo. 


Proof. Notice that md{o') is asymptotically equivalent to 

[2(1 + e^'") + hcT{l - hb^d^ 

as d —>■ oo, the roots of which are given by (l2^ . □ 


Notice that dsup is independent of the system parameters: it depends on h 
only. In our example we have h = 0.1, which gives dsup = 23.99 and agrees with 
Fig. [5l 

The notion of impedance matching^ from the theory of transmission lines, 
suggests the choice d = /cp in the scattering transformation (ED. The ratio¬ 
nale behi nd this choice is that, in a real transmission line, it avoids wave re¬ 
flections Niemever and Slotind . 12004 1. From a frequency-domain perspective, 
the immediate advantage of this choice is that the characteristic quasipolyno¬ 
mial m simplifies substantially, as the coefficients become 


P2{s) = 2kp , 

pi(s) = 2kp{a + bkp) -|- (1 - e~^’')ki , 

Pq{s) = (1 -I- e~^'‘)bkpki -|- (1 - e~^‘')aki . 

The principal coefficient becomes constant, i.e., the closed-loop system is of 
retarded nature instead of neutral, which obviates the need to verify the stability 
of the difference operator given by the additional constraint (ED- 
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The restriction d = kp corresponds to the plane shown in Fig. |4(a)[ Notice 
that, while it is a reasonable choice because of the arguments given above, the 
plane fails to intersect many u-stability surfaces. Thus, the choice d = kp is not 
optimal in the sense that it obstructs the achievement of cTsup-stability. However, 
notice from Fig. [S] that the normalized minimal gain, kp{a^)/d, converges to a 
constant value as d —> oo. This fact suggests the more general linear relation 


d — ; 


(23) 


where C > 0 is a new design parameter. In the reminder of this section, we 
will show that there exists a privileged value of <(, which we call ((min and that 
optimizes the bound on tr. 

Let us first compute the (T-stability boundaries on the plane given by (USD- 
The boundaries associated to the real roots s = —a are simply given by (USD 
and (PSI) . On the other hand, solving (ITS)) with d = (^kp is slightly more difficult, 
since A(lo) and B(lo) now depend on kp. It follows from Cramer’s rule that 

\A{uj)\kp + d\{A 2 {uj) i?(a;))|= 0 , ( 24 ) 

where ^2(0;) stands for the second column of A{oj) and | • | stands for the 
determinant. Substituting (ESI) into (ESI) gives the equation 

|H(a;)| + (( I (^2(0;) B[u})^ | = 0 . 

Developing this equation explicitly leads to the second order polynomial equa¬ 
tion C2b‘^C^kp + cibCkp -I- Co = 0 with 

C2 = ((C “ l)e^^°’ — (C -I- 1 ) — 26^“^ cos{huj)) oj + si'n{huj)a , 

Cl = 2 ((—(C ~ “ (C + 1 )) + 2e^°’ {C,acos{huj) +ojsm{ha))) tj , 

Co = -I- a^) ((e^^'’"(C — 1) — (C + 1)) w + 26^^“^ {ujcos{huj) — ((jsin{huj)) ) . 


The roots of this polynomial can be computed explicitly, 
computed as 


h = ^12^(w) {dBi{oj) - Aii{uj)kp) 


Finally, ki can be 


with 

Fig. 6(a) shows the cr-stability regions for the plane d = kp. Notice that, 
as in Fig. 2(a)[ the boundaries are not given by closed curves. This is a clear 
advantage with respect to the case in which d is constant (cf. Fig. |4(b)[ ). Un¬ 
fortunately, the choice d = (^kp still imposes a bound on the exponential decay, 
as it will be shown shortly. 


Definition 1. Let 77sup be the positive solution of 


2(1 -h e''"“p) -h r?sup(l - e''"“p) = 0 


(25) 


and let ((niin be defined as 


(I -I- 

2??supe''™p — (1 -I- e’'="p)(I — e'^^pp) 
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(a) With C = 1- The exponential decay is bounded by = 12.78. 
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Figure 6: u-stability regions on the plane d = Qkp. Communication channel with 
delay, h = 0.1, and scattering transformation. The diamond markers correspond 
to the minimal gains that ensure a given cr. 





16 














Proposition 4. Consider a plant (HD in closed-loop with a PI controller m 
satisfying kp > 0 and ki > 0. A communication channel with round-trip delay 
h > 0 stands between the system and the controller. Set d = C,kp with C > 0 
apply the scattering transformations m at the channel end points. 

(i) The least upper bound on the exponential decay is 




if 0<(<(. 


cr > 0 such that m({ha) =0 if Cmin < C 


where 

m^{ri) = (1 + e^) (C(l - e'') + 1 + e'') - 2C,r]e^ . 


(26) 


(ii) The minimal kp assigning a given > a > aj^ is a root of the second 
order polynomial C 2 h^Cf^kp + cibfkp + cq with 

C2 = (1 + e'^'^) (C(l - + 1 + e'“") - 2 (:hae^^ , 

Cl = 2(1 + e'^'^) (C(l + + 1 - a, - 4 Cae'‘" , 

Co = ((1 - (C(l + + 1 - e'*") + 2 Chae'^^) al , 

while the minimal ki is given by 

_ (1 + e^'^)(((&fcp + Co-) + (1 — e^'^)(^^6A:p + Co-) , 

* “ (1 + e^'")C6fcp + (1 - e^'^)aa- ^ ^ ' 


Remark 3. Unlike the bounds given in Propositions [T] and [2l the bound 
is not achievable in practice, as it requires infinite gains. It can, however, be 
approached arbitrarily close. 

Proof. The coefficients of (fT4l) take the form 

P2{s) = (1 + e“'*®)Cfcp + (1 - , 

pi(s) = (1 + e~^’'){bkp + a)Ckp + (1 - e~^'‘){bC,'^kp + akp + h) , 

Po{s) = (1 + e~^'')bhCkp + (1 - . 


The formulas for computing the minimal gains are obtained by setting (1141) and 
its derivative equal to zero. This proves (ii). 

The maximal exponential decay approaches its maximal value, cr^((C), as 
kp and ki tend to infinity, that is, as C 2 approaches 0. Notice that C 2 = 0 is 
equivalent to mQ{ha) = 0. For hxed there exists an 77 > 0 such that m(^{rii) — 0 
if, and only if, 

Cmin ^ C ■ 

To see this, set m(^{r]) = 0 and write ^ as a function of 77 : 


(1 + e^)^ 

277e')-(l + e'')(l-e’') ' 


(27) 
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To find the lower bound we solve mQ{r]) = 0 and 

^771^(77) = - 2 e ^ [C(e’' + ?7 + 1)) + 1 + o'*] = 0 

simultaneously for ry. This gives the implicit equation (1251) . The lower bound 
on C is finally found by substituting rjsn-p in (l27)l . that is, it is given by Cmin- 
Thus, mc^{ha) = 0 is solvable only if ^min < C- For 0 < C < Cmin, condi¬ 
tion dm) determines the maximal exponential decay. When d = Ckp, condi¬ 
tion El) reads 


This proves (i). □ 

The following theorem is important from a practical point of view, since it 
gives an objective choice for the free parameter in the scattering transformation. 

Theorem 1. Consider a plant ([T]) in closed-loop with a PI controller ([2|) sat¬ 
isfying kp > 0 and ki > 0. A communication channel with round-trip delay 
h > 0 stands between the system and the controller. Set d = Cmin^p o-nd apply 
the scattering transformations ED at the channel end points. The least upper 
bound on a, is equal to CTsup- 

Proof. To compute Cmin we first solve (EH), which is equivalent to (1^^ with 
g = hasnp- □ 

Remark 4. The theoretical limit on the exponential decay is independent of the 
plant parameters. It depends only on the delay h and it is given by 

_ Vsnp _ 2.3994 
h ~ h ■ 

The optimal parameter Cmin = 0.8336, for which a can be set arbitrarily close 
to CTsup, is independent of the plant parameters and the delay. 

In order to assist the reader, the necessary stability condition of the difference 
operator is plotted in Fig. [7] with d = C,kp along with the solutions of equation 
= 0 in (l26l) . As it results from the above analysis, C is indeed bounded 
from below by Cmm- Following this observation, the stability of the difference 
operator is always satisfied in the interval Psup/h > cr > 0 and therefore, the 
(T-stability of the overall control-loop is finally established. 

The (T-stability regions on the plane d = Cmin^p are shown in Fig. |6(b)| 
(compare with Fig. |6(a)] ). Notice that is almost twice as large as 

the value obtained using the recipe usually found in the literature. The plane 
is also shown in Fig. |4(a)| 

In view of the preceding remarks we propose the following design procedure: 

(i) Set d = Cminfcp- This choice is universal, in the sense that it does not 
depend explicitly on the plant parameters nor the delay. 


ha < In 


1 + C 


1-C 
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Figure 7: Pairs (r/,() solution of mdr]) = 0, 77 > 0 (in blue). Bounds given 
by uni) (in green and red). Thick lines correspond to the least upper bound on 
the exponential decay, multiplied by h. 


(ii) Choose a desired cr such that rjsup/h > a > 0. This choice requires 
knowledge of the round-trip delay only. The choice of a should take into 
account the limits imposed by the actuator, i.e., large a may generate 
saturation. 

(iii) Use the corresponding minimal gains, as in Proposition 31 This choice 
requires knowledge of the round-trip delay and the plant parameters. 

4 Concluding Remarks 

A simple instance of the use of the scattering transformation in passivity-based 
control has been analyzed from the classical (i.e., frequency-domain) perspec¬ 
tive of time-delay systems. The exponential decay rate of the closed-loop system 
was chosen as a criterion to asses the performance of different control schemes. 
This criterion leads to an optimal choice on the design parameter of the scat¬ 
tering transformation. Quite remarkably, the optimal choice is independent of 
the plant parameters and the delays in the communication channel. With the 
optimal choice, it is possible to attain exponential decay ratios that are almost 
twice as large as those obtained by setting the design parameter at the value 
suggested in the literature. 

A theoretical limit on the exponential decay rate has been found. The limit 
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does depend on the total delay but, quite remarkably as well, it is independent 
of the plant parameters. 

Distinct qualitative features, such as the shape and extension of the a- 
stability regions or the neutral and retarded nature, have been identified for 
different scenarios: a closed-loop system without delays, with delays, with and 
without scattering transformation. This furthers our insight on the effect of 
delays and the scattering transformation used to remedy it. 
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